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Abstract 

We use the relativistic light-front quark model to show that both strong and 
radiative D* decays are in good agreement with the 1992 CLEO II results. In 
particular the coupling for D* — > Dir is consistent with the experimental upper 
limit. The key point is the relativistic treatment of the quark spin. 
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1 Introduction 



The D* decays has been under both theoretical and experimental investigation for 
quite a long time ||, g, g, |, |, § Q |, fj. The strong D* decays, D* -> JJtt, can be 
described by the lowest-order (in external momentum) effective Lagangian 

C = (^)DV ■ ^TT^ , (1) 

where Z) and Z)* are the D and D* isodoublet fields, and r is the isospin matrices. 
The decay width of D* + — > D°tt + , for example, reads 



_ £ 2 |P^| 3 



(2) 



The radiative 7J* decays 7J* — > /J7, on the other hand, can be described by the 
decay amplitude 



M(D* - 7J 7 ) = e/i e aiATA C^P^Pd , (3) 

where e/i/2 is the transition magnetic moment, and £ a and e u are the polarization 
vectors of the photon and the D* meson, respectively. The radiative decay width is 



<«*> w • (4) 



In many theoretical studies 0, [5| people use the non-relativistic quark model to 
calculate the radiative D* decays while they use various modifications of the 577(4) 
symmetry method for the strong D* decays. Obviously, the treatment of the strong 
and the radiative D* decays is not on the same footing. 

However one can also use the quark model to calculate the strong decays. Through 
the PCAC relation one can relate the strong coupling g to the axial-current form 
factor A (0) 

g = A° D *(0). (5) 
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Here Aq is one of the four form factors used to describe the matrix element of a vector 
meson V decaying to a pseudo-scalar meson X by the quark transition Q' — > Q: 



; le^e p x p^-2m v - 

mx + my q 



(X{p x )\ Q 7/1 (l - 75 )Q' |y(pv,e)) =f^rV«^ W P^ -2m v ^ J A V A)(g 2 

(6) 



(m x + m v )e^A x { q 2 )- l " ' P * J (p* + PvTA 2 {q 2 ) - 2m y ^|^g^ 3 (g 2 ) 

mx + vny q z 



where e is the polarization vector of the vector meson, q = Pv — Px, and 

2 mx + my 2 my - m X . , 2 \ a / r\\ a / r\\ f-\ 

Mf) = -j^-Mq 2 ) ^— Mq 2 ) , MO) = MO) . (7) 

In Eq. (|^), Aq D *(0) is the A (q 2 = 0) form factor for D* — > D induced by the quark 
transition u — > d. In the non-relativistic quark model, one obtains g — 1 [|Kj. In a 
rather different model, the BSW model JTTJ, one has 



3 = < D *(0) = Jdx d 2 k ± 0^(x,k ± )0 D »(x,k ± ) , (8) 

where 4>(x, kj_) is the momentum wavefunction. Since the wavefunctions of the D 
and .D* are quite similar, this model also has g — 1. In fact, in the heavy c-quark 
limit 4>d*(x, kj_) and 4>d{%, kj_) should be exactly the same, due to the heavy quark 
symmetry, and hence in the BSW model, g = 1 in this limit. These results of g are, 
however, larger than the upper limit of g obtained from the ACCMOR measurement 
IV+ < 131keV §: 

9 < 0.7 . (9) 

There have been some recent studies on the D* decays || [7|, §, |J. In ref. || 
the radiative D* decays are treated with the non-relativistic quark model but the 
strong coupling g is related to the quark-level axial-current form factor g u ^ = ^ ; this 
is obtained from fitting the measured nucleon axial-current form factor g° ucleon = 
1.25 = ^g^. Thus the coupling g is g = g u ^ = ^. In refs. |7|, H heavy-meson 
chiral-perturbation theory [|T2 , |10[1 is used and all the D* decays can be described 



by two parameters. These parameters are the couplings of the chiral Lagrangian 



and can therefore be determined only by fitting the experimental data. In a recent 
study |J , the radiative decays are calculated using the vector-dominance hypothesis 
while the strong coupling g is extracted from the data on the decay D —>■ Treu e using 
heavy-meson chiral-symmetry relations. The central value obtained for g is g ~ 0.4. 
In this study the uncertainties associated with assumptions on the /+ form factor of 
D — > 7rez/ e and on the decay constant fjj, etc., could be large. A value of g ~ 0.4 is 



also obtained |T3[ in a numerical solution to another type of relativistic quark model 
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The failure of the non-relativistic quark model and the BSW model in explaining 
g < 0.7 has its roots in the non-relativistic treatment of the quark spin in these 
models. In this letter, we present a calculation of the strong and radiative decays of 
D* using the relativistic light-front quark model [15. [TjJ 17]. In our calculation we 



are able to treat on the same basis both the strong and radiative decays. Due to the 
relativistic treatment one can see unambiguously how the strong coupling g can be 
consistent with the limit in Eq. (|9|). Our calculated branching ratios for the strong 
and radiative decays are in good agreement with the 1992 CLEO II measurement 



2 The light— front relativistic quark model 

The relativistic light-front quark model was developed quite a long time ago and 



there have been many applications [l_g, [16, [T7[ • Here we give a brief introduction 



A ground-state meson V(Qq) with spin J on the light front can be described by 
the state vector 



\V(P, J 3 , J)) = J d 3 Pl d 3 p 2 8(P - Pl - P2 ; 



£ ^ j3 (P,p 1 ,p 2 ,A 1 ,A 2 )|Q(A 1 , Pl ) g(A 2 ,p a )> , (10) 

Ai,A2 
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In the light-front convention, the quark coordinates are given by 

p+ = Xl P + , -pt = X 2P + , Xi + x 2 = 1 , < Xi j2 < 1 , 

Pl_L = XiP± + k ± , p 2± = x 2 P± - k ± . (11) 

The quantities x\ >2 and k^ are invariant under the kinematic Lorentz transforma- 
tions. Rotational invariance of the wave function for states with spin J and zero 
orbital angular momentum requires the wave function to have the form Jl^, |TJJ] (with 

X = X\) 

tt J - Ja (P, Pl , P2 , A 1; A 2 ) = R J ^(k ± , A 1; A 2 )0(x, ki), (12) 
where <p(x, kj_) is even in kj_ and 

R J > J *(k ± ,\ 1 ,\ 2 ) =^(A 1 |4 / (k ± ,m Q )|A)(A 2 |4 f (-k ± ,m 9 -)|A')C J ' j3 (i A;i A').(13) 

A, A' 1 1 

1 1 

In Eq. (jl3|) , C J ' J3 (2, A; A') is the Clebsh-Gordan coefficient and the rotation 
i?ju(k_i_, rrii) on the quark spins is the Melosh rotation [IS|: 



T> n m + x { M -itr -(n x k ± ) 

.flAf(ki, mi) = — , (14) 

(m, + x.Mq) 2 + ki 



where n = (0, 0, 1), a is the Pauli spin matrix, and 



„ m 2 n + k^ m? + k^ 
M n 2 = — S ^ + —2 i . 15 

Xi x 2 



The spin wave function R J,j3 (k±, Ai, A 2 ) in Eq. (13) can also be written as 



i^ j3 (k ± , A 1; A 2 ) = X \Ai^ m Q) S ' V3 4T(-k±,m 9 -) XA2 



where S ,J,J3 is defined by 



5 J - J » = j:iA><V|C jr ' J -(i A;i A'). (17) 



A, A' 



For the pseudoscalar and vector mesons 

i<J 2 

71 



S 1 - = ^= . 
\/2 



(18) 



The explicit expressions of [7 1 " 73 in Eq. (pi) can be found in ref. |TT| 



The matrix element of a vector meson V(Q'q) decaying to a pseudo-scalar meson 
X(Qq) is 



(X(p x )\QTQ'\V(p v ,J 3 )) = Jdx d 2 k ± 



dx d k i 



Ai,A 2 



Tr 



x 



(19) 



where is defined by 



U Q F M Q' = xl U r Xj 



(20) 



In Eq. ([19]), we choose g + = so q 2 = — c[ 2 ± . In contrast to Eq. (|T9|), the matrix 
element in the BSW model is given by 



(X(p x )\QTQ'\V(p v , J 3 ))= / d 2 k x ^^Tr IS^UrS 1 ^ 

J x 1 



(21) 



Eq. ( |T9"D is in fact expected to be valid only for "good" currents such as T - 
7 + , 7 + 75, .... There are contributions other than the one given in Eq. ([II]) if 



current is not a "good" current [[16 



Using the "good" currents T = 7 + 7s,7 + and Eq. (|T9"D, we obtain the following 
expressions for the form factors A (0) and V(0) of the Q' — > Q transition 



A o (0) 
V(0) 



dx d 2 k±4>* x (j)v 
(^+ki)(^+ki 

<ix <i 2 k_|_03c0y 
(^+ki)(^+ki 



A v ^ + (2x - l)ki + 2(m y +mq)(l-x)k3 

W v 

--(m x +m v )(l-x) 



A x +-± + (1 - s)(mg - mgOki^y 



(22) 



(23) 
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where 



Ay = x rrig + (1 — x)mqi , A x = x rriq + (1 — x)mQ 

W v = Ml + m Q + m- q , 9 V = (^f)/4>v • (24) 

Here, corresponds to Eq. (|1|) for the meson V. 

The meson wave functions (p(x, kj_) are model dependent and difficult to obtain; 
often simple forms are assumed for them. One possibility is to use the wavefunction 
adopted in [Tl| 

k 2 



x, kj_) = NJx(l — x)exp 



M 2 
2w 2 



1 itlq — m?l J '\ ' ri> 



2 2M 



2 



2w< 



7rur 



(25) 



where M is the mass of the meson. In ||TB| a Gaussian type of wavefunction was 
used, 
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where k z is defined by X\M Q = Eq + k z with Eq = y k^ + fc 2 + tuq. The parameter 
lj in both Eqs. (|25| ) and (|26|) should be of the order of Aqcd- As we will see, the 
results are not too sensitive to the choice of wavefunction. 

It is interesting to look at the general behavior of the wavefunction (p(x, kj_) in 
the limit of tuq — > oo. The distribution amplitude / rf 2 k^0(x, kjj of a heavy meson, 
as is well known, should have a peak near x ~ xq = As vtiq becomes larger the 
width of the peak decreases and xq comes closer to 1. The wavefunction 0(x, kj_) 
vanishes if k^ > > Aq CD and in the tuq — > oo limit peaks at Xq — ^ 1 cLS the width of 
<f)(x, kj_) — > 0. Both wavefunctions listed have this feature. 

3 The D* decays 

We now calculate the D* decay rates in the light-front quark model. For D* — > Dtt, 
we first look at the form factor Aq(0) of the transition Q — > Q. In this case Eq. (p2]) 
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reduces to 



where 



A$ Q (0) = Jdx d 2 k ± ^0 y T Ao (27) 



A 2 + (2s - l)ki + AmQil w ~ x)K 
Ta = ^-^ v - (28) 



and 



A = x rriq + (1 — x) tuq . (29) 

Generally Ta < 1 since it comes directly from the Melosh rotation. (For no 
Melosh rotation one would get Ta — 1.) How far Ta deviates from 1 depends on 
the relation among tuq, rriq and kj_, where the average value of is determined 
by the scale parameter in the wavefunction such as the u in Eqs. (P5|) and (|26|). 
Usually as < k^_ > becomes smaller Ta becomes closer to 1. If k^ = in Eq. (fffi) 
then Ta = 1. 

Let us consider two limiting cases for the transition quark Q. If the transition 
quark Q is infinitely heavy while the spectator quark is light, x — > 1, and therefore 
the third term in the numerator of Ta vanishes. Thus 

T Ao -> 1 and A^ Q (0) -> / dx d 2 \t L <p* x <p v -> 1 , (30) 

as required by the heavy quark limit. Note in this case q 2 = q 2 max = 0. If, as in the 
decays D* — > Z>7r, the transition quark Q is light while the spectator quark is heavy, 
then in Eqs. ( p7|) and (p8[) , x — > as m c — > oo. In this case the third term in the 
numerator of Ta also vanishes, but now 

A 2 - k 2 

t Ao - -^rp^- . (31) 

Since both k^ and A 2 are of the same order, one can expect considerable deviation 
from 1 for A (0) in this case. 



Obviously, the situation for g = Aq D * (0) is quite close to the second case dis- 
cussed above. Thus, one can clearly see why the light-front model gives a deviation 
of g from 1. 

To calculate the radiative decays D* — > Dj, one can express \i in Eq. in 
terms of the form factor V(0) 

/i = [e c ^ c (0) + e,V*(0)]( I ) (32) 

where q = u or d and e c = ^, e u = ^, = — tj. V^ Q (0) is the V(0) form factor 
of the transition Q —+ Q (Q = c,u or d). V^^O) in the light-front model can be 
obtained from Eq. (|23|) 

V*(0) = (m p + mo .)J " X d ^°- (1-x)[a + ^_] . (33) 

To calculate V cc (0), V^ q (0), and g, we use the two wavefunctions given in Eqs. 
(25) and (EH) as the wavefunctions for the Z) and D* mesons. We assume m„ = 



rrid = m and expect that the light-quark masses m in the wavefunctions to be 
around 0.25 ~ 0.30GeV. Similarly the scale parameter u in both wavefunctions 



should be around 0.40 ~ 0.50GeV. These values are similar to those taken in |TT 
and . 



The dependence of V^^O) on the quark mass rriQ is similar to that in the non- 
relativistic quark model 0, [20]: V QQ = { mD \ mD * )^/jf^ (jh) ■ That is V Qq (0) 
decreases as rriQ increases and if vtlq — > oo, V^(0) — ► 1, as required by the heavy 
quark relation. 

In the figure we show the dependence of the coupling g on the light quark mass 
m and u in the wavefunctions Eqs. (p5| ) and fl26"|) . We choose m c = 1.6GeV; the 
coupling g has little dependence on m c . One can see that for a fixed lu, g increases 
as m increases. For m between 0.25 ~ 0.30GeV and u between 0.40 ~ 0.45GeV, g is 
between 0.55 ~ 0.65. To avoid adjusting parameters, we choose m = 0.25GeV, and 
uj = 0.4GeV for the calculation of decay rates. There is good agreement between the 
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calculated branching ratios and experiment. With this set of parameters g = 0.6 
using Eq. (|25|) or Eq. (|26|) . In the table , we show the corresponding results for 
D* — > Dtt and D* — > D7. We also show how the branching ratios change with a 
different choice of mass. Finally, we give the results of our calculation for the decay 
D* + — > (for m s = 0.40GeV and m s = 0.50GeV). This decay is usually assumed 
to be the dominant mode. In a recent paper |^TJ it has been argued that through 
isospin symmetry breaking there can also be a significant correlation between the 
branching ratios Br(D* — > D s tt°) and Br(D* + —> D + j). 

4 Conclusion 

In this letter we have presented a calculation of both strong and radiative D* decays 
using the relativistic light-front quark model. Our results are in good agreement 
with the 1992 CLEO II measurement. In particular the strong coupling g is con- 
sistent with the experimental upper limit Td*+ < 131keV. The key point is the 
relativistic treatment of the quark spin. The fact that the relativistic treatment 
is essential in explaining g < 1 is reminiscent of the similar situation in the nu- 
cleon axial-current form factor g^ cleon where the non-relativistic quark model gives 
(?nucieon = -| and a relativistic treatment can get g^ ucleon naturally down to the mea- 
sured value . 
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Figure: 

The dependence of the coupling g on the light-quark mass m and to in Eqs. (f25|) 
and (pr^). The solid lines are the results of Eq. ( p6|) and the dashed lines use Eq. 
(|25|). The upper solid and dashed lines correspond to u = 0.40 GeV, and the lower 
solid and dashed lines are for uj = 0.45 GeV. 
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Table. Numerical results for D* decays. Also listed are the experimental results of 
CLEO II B and PDG §. All branching ratios are in %. 



Decay 




Br «) 


g r a) 


Br fe ) 


Br 6 ) 


Br 


Br 






Eq.(H) 


Eq.(|26|) 


Eq-dH) 


Eq.© 


CLEO II 


PDG 


D*+ - 


>D°7I + 


68.2 


68.0 


68.1 


67.7 


68. 1± 1.0± 1.3 


55± 4 


D* + - 


► D + TT° 


31.6 


31.5 


31.6 


31.4 


30.8± 0.4± 0.8 


27.2± 2.5 


D* + - 


► D + 1 


0.3 


0.5 


0.4 


0.9 


1.1± 1.4± 1.6 


18± 4 


D*° -> 




75.2 


73.0 


70.5 


66.2 


63.6± 2.3± 3.3 


55± 6 


D*° -> 




24.8 


27.1 


29.5 


33.8 


36.4± 2.3± 3.3 


45± 6 






r(keV) a ) 


r(keV) a ^ 


r(keV) 6 ) 


r(keV) b > 










Bq.® 


Eq.(|2§ 


Eq.®) 


Eq.© 






D* + - 


> total 


122.8 


123.0 


104.9 


102.1 






L>*° -» 


total 


74.0 


76.3 


67.4 


69.5 






D*+- 


► D+7 


0.1 


0.1 


0.2 


0.3 






a): uj = 


0.4GeV, 


m = m u 


= m d = 0.30GeV, m s = 


0.50GeV. 






b): u = 


0.4GeV, 


m = m u 


= m d = 0.25GeV, m s = 


0.40GeV. 
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